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> : ABSTRACT 

cn 

Q\ • We have developed a numerical model for the temporal evolution of parti- 

^ ' cle and photon spectra resulting from nonthermal processes at the shock fronts 

O ■ formed in merging clusters of galaxies. Fermi acceleration is approximated by 

^ ■ injecting power-law distributions of particles during a merger event, subject to 

constraints on maximum particle energies. We consider synchrotron, bremsstrah- 
2 ■ lung, Compton, and Coulomb processes for the electrons, nuclear, photomeson, 

^ ■ and Coulomb processes for the protons, and knock-on electron production during 

the merging process. The broadband radio through 7-ray emission radiated by 
nonthermal protons and primary and secondary electrons is calculated both dur- 
ing and after the merger event. Using ROSAT observations to establish typical 
parameters for the matter density profile of clusters of galaxies, we find that typ- 
ical merger shocks are weak and accelerate particles with relatively soft spectra. 
We consider the prospects for detecting nonthermal radio and 7-ray emission 
from clusters of galaxies and implications for the origin of ultra-high energy cos- 
mic rays and the diffuse 7-ray background. Our results suggest that only a few 
of the isotropically-distributed unidentified EGRET sources are due to shocks 
formed in cluster mergers, and that only a minor contribution to the diffuse ex- 
tragalactic 7-ray background can originate from cluster merger shocks. Cluster 
merger shocks can accelerate protons to < 10^^ eV for the standard parame- 
ters considered here. We predict that CLAST will detect several cluster mergers 
and, depending on the mean magnetic fields in the intracluster medium, the Low 
Frequency Array could detect anywhere from several to several hundred. 
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1. Introduction 

According to the hierarchical merging scenario, cold dark matter halos evolve to form 
larger structures by merging with adjacent dark matter halos. Within dark matter halos, 
baryonic matter condenses to form clusters of galaxies. A cluster merger event results from 
the interaction of galaxy clusters during the merger of cold dark matter halos. The grav- 
itational potential energy available in a cluster merger event involving halos with masses 
~ lO^^M© is ~10^^ — 10^^ ergs. As two clusters of galaxies merge, the infall velocities can 
exceed the sound speed of the ICM. As a result, a shock front will form at the interaction 
boundary between the clusters. First-order Fermi acceleration at the shock front produces 
a population of nonthermal, relativistic particles. Relativistic electrons are detected from 
their synchrotron radio emission or from 7-rays due to Compton-scattered cosmic microwave 
radiation. Nonthermal protons are detected through 7 rays emitted from secondaries formed 
in nuclear production processes, including the 7r° decay signature at 70 MeV. 

Rich chisters contain a hot and tenuous ionized intracluster medium (ICM), with ob- 
served temperatures Tx ~ 5 — 10 keV, sound speeds ~1000 km s""*^, and thermal bremsstrah- 
lung luminosities Lx ~ 10^^ ergs s~^ between 2 and 10 keV. In addition to the luminous 
thermal component present in these clusters, there is a growing body of evidence supporting 
the presence of nonthermal distributions of particles in cluster mergers (Eilek 1999; Feretti 
et al. 2000). Deep radio observations of clusters of galaxies indicate the presence of extended 
diffuse emission not easily associated with an optical counterpart. These diffuse radio fea- 
tures are commonly classified as either radio halos, or radio relics (which are also called 
periphery halos). Radio halos mimic the observed X-ray profiles and are characterized by 
their central location in the cluster and by a highly disorganized magnetic field, and are 
generally thought to be a consequence of a merger event (Feretti et al. 2000). Radio relics, 
usually found on the periphery of the cluster, are characterized by highly organized mag- 
netic fields and often display filamentary structures. We focus here on radio relics, which 
are thought to result from synchrotron emission emitted by electrons directly accelerated 
at shock fronts. By contrast, radio halos might be due to reacceleration of rehc electrons 
by magnetic turbulence or enhanced magnetic turbulence arising from the cluster merger or 
motions of the galaxies (Brunetti et al. 2001; Ohno, Takizawa, & Shibata 2002). 

The first radio halo in a cluster of galaxies was detected from the Coma Cluster (Large 
et al. 1959). Many halos and rehcs have since been found in a number of other clusters, 
including A754 (Kassim et al. 2001), A2256 (Berrington, Lugger, & Cohn 2002), and others 
(see Govoni et al. (2001); Slee et al. (2001) for recent observations). Typical radio powers 
from the radio halos and relics are at the level of lO'^'^-lO*^^ ergs s^^ (Giovannini & Feretti 
2000). These clusters show a favorable correlation (Kassim et al. 2001; Berrington, Lugger, 



-3- 



& Cohn 2002) between the existence of diffuse radio features and recent or on-going merger 
activity. Detection of these features is largely due to sensitivity limitations; the number of 
detected diffuse radio features has increased with the improvement in sensitivity of radio 
telescopes. 

The synchrotron evidence for nonthermal electrons in clusters of galaxies indicates that 
an acceleration mechanism must exist in the environment of the host cluster with sufficient 
power to produce the observed emission. Several mechanisms have been proposed to explain 
the correlation of the radio halo and relic features and the recent or on-going merger activity. 
These mechanisms often require the presence of shock waves or magnetic turbulence to 
accelerate particles via the Fermi acceleration process (Schlickeiser, Sievers, & Thiemann 
1987; Tribble 1993; Kang, Rachen, & Biermann 1997; EnBhn et al. 1998; Blasi 2000; Miniati 
et al. 2001). Other theories include adiabatic compression of fossil radio plasma by a cluster- 
merger shock wave (EnBlin & Gopal-Krishna 2000). 

Optical surveys show that approximately 30-40% of clusters of galaxies display evidence 
for the presence of substructure (Forman et al. 1981; Geller & Beers 1982, and others). This 
internal structure is often interpreted as a subset of galaxies merging with a larger cluster 
of galaxies. These internal structures indicate velocity differences near or greater than the 
expected sound speed of the intracluster medium (ICM). Observed velocities typically range 
from ~ 1000 to 3000 km s~^ which are consistent with values expected from parabolic orbits 
(Oegerle & Hill 1994). 

Emission from nonthermal particles will also appear at EUV and hard X-ray energies as 
a power- law excess. While the EUV emission seen in the Coma cluster (Lieu et al. 1999), and 
possibly also A2199 and A1795 (Lieu, Bonamcntc, & Mittaz 1999) may have a cool thermal 
origin, it is unlikely because of the extreme mass requirements of cool gas. It is more likely 
that the EUV emission has a nonthermal origin (Hwang 1997; Enfilin & Biermann 1998; 
Mittaz, Lieu, & Lockman 1998; Blasi & Colafrancesco 1999; Atoyan & Volk 2000). Excess 
hard X-ray (HXR) emission has been reported from the clusters A1656 (Fusco-Femiano et al. 
1999; Rephaeh, Gruber, & Blanco 1999), A2256 (Fusco-Femiano et al. 2000), A3667 (Fusco- 
Femiano et al. 2001) and possibly A2199 (Kaastra et al. 1999). Estimated luminosity of 
the HXR emission are fslO^^ ergs s~^. Thermal origins of this HXR excess emission would 
require unrealistic temperatures greater than 40 keV, and so is thought to be caused by the 
presence of relativistic nonthermal electrons (Fusco-Femiano et al. 1999). 

Numerical models of merging clusters of galaxies (Ricker 1998; Takizawa 2000; Ricker 
& Sarazin 2001; Miniati et al. 2001, and others) have treated the development of shocks as 
a result of the cluster merging process. Given the presence of thermal ionized particles in 
the vicinity of these shocks and a cluster magnetic field, electrons and ions will be acceler- 
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ated via the first-order Fermi process (Bell 1987; Blandford & Ostriker 1987). Numerous 

attempts have been made to model the emissions from nonthermal particles (Colafrancesco 
& Blasi 1998; Fujita & Sarazin 2001; Petrosian 2001; Miniati et al. 2001) produced by these 
cluster merger shocks, but this study differs from previous attempts in that we accurately 
model the diffusion of particles in energy space due to Coulomb interactions, and allow for 
a variable injection rate that depends on the environment local to the shock front. Recent 
studies by Liang, Dogiel, & Birkinshaw (2002) have addressed the balance of Coulomb losses 
and diffusion in energy space, but do not address the injection and diffusion of particles in 
energy space along with a variable source function of nonthermal particles. In our treatment, 
we follow particle energies up to 10^^ eV, and calculate the bremsstrahlung, Compton, syn- 
chrotron, and TTo 7-rays from p-p collisions. Because the shock front lifetime is a significant 
fraction of the age of the universe, we also follow the changing environment due, for example, 
to the changing CMB energy density. In addition we accurately model nonthermal electrons 
and protons up to ~10^^ eV, though we find that limitations on particle acceleration make 
it difficult to produce protons above ~ 10^^ eV in cluster merger shocks. From the particle 
distributions, we calculate nonthermal photon spectra for energies up to ~10^ GeV. We also 
include the effects of secondary production on the nonthermal photon spectra. 

The physical processes and the temporal evolution of the particles are described in 
section 2. The results of the simulations are presented in section 3. A comparison of the 
photon spectra with observed clusters, and the potential of detecting these shocks with 
space-based sateUite observatories is discussed in section 4. 

2. Models 

We have developed a code to calculate the time-dependent particle distribution functions 
evolving through radiative losses for electrons and protons accelerated by the first-order 
Fermi process at the cluster merger shock. The code was originally adapted from a supernova 
remnant code (Sturner, et al. 1997), and apphed to the specific case of the cluster merger 
scenario. We will break up the discussion of the problem in to four subsections: cluster 
merger dynamics, nonthermal particle production, temporal evolution, and the production 
of the photons. 
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2.1. Cluster Merger Dynamics 

2.1.1. Evolution of the Merger Event 

Two clusters are assumed to be the dominant gravitating masses in the local region, and 
variations of the gravitational field due to surrounding masses are assumed to be negligible. 
A cluster accretes another cluster which falls from a distance d. This distance can be 
approximated by the maximum separation of two point masses which move apart beginning 
at the moment of the big bang, follow an elliptical orbit, and merge at time t„j. Given 
this assumption, a cluster of mass Mi = IO^^A'/isMq accretes a subcluster of mass M2 = 
lO^^misMo at time t„i as approximated by Kepler's third law. Consequently, a cluster that 
merges at a time t^ fell from a maximum separation approximated by 



d 



2G{Mi+M2)tl 
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3 
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(Ricker & Sarazin 2001). This equation allows us to calculate the magnitude of the energy 
pool available for particle acceleration. The total orbital energy £ of the merging system, 
assuming a zero impact parameter, is given by 

~ 2l/3(Ml+M2)V3 ytrr^) 

/ \ -2/3 

^ ^ -1-9 >< 10'' (mSS^ (10%^) ergs. 

We use the terminology that orbits are bound if the total energy £ < 0; otherwise, they are 
unbound. It is customary to use the Keplerian analogy of elliptical to represent orbits where 
£^ < 0, and parabolic for orbits where £ = 

To approximate the merger velocity of the system, we relax the assumption that the 
accreting cluster is a point mass. Instead, we assume the ICM number density of the accreting 
cluster follows the spherically-symmetric, isothermal beta model given by 



2i/3(Mi+M2)V3 yt„,) 



nicuir) = riQ 
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(3) 



The quantity no is the central number density of the ICM, and Vc is the core radius, which 
characterizes the spatial scale over which the density changes from a roughly constant value 
at r < Tc to a power- law density dependence at r > Tc. The gravitational dynamics are 
dominated by the dark matter distribution of the more massive cluster. Assuming that the 
dark matter distribution follows a similar profile as the ICM, and that the merging cluster 
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is a point mass, we calculate the cluster merger velocity and radial separation by solving the 
differential equation 

dt'^ dri2 

The radius ri2 is the separation of the center of masses of the two clusters, Mi(ri2) is the 
mass interior to the subcluster at radius ri2, M2 is the mass of the merging subcluster, and 
= M^^M reduced mass. The gravitational potential $(Mi,ri2) is defined by 



4,(M,r)=-^. (5) 



The mass interior is defined as 



M{r) =A7rU dx x^mcuix) (6) 
Jo 

where ^rn is a normalization constant chosen so that Mi(ri2) is equivalent to the total mass 
of the dominant cluster at its maximum radius Ri. Obviously, the potential satisfies the 
requirement $(M, r — > 00) = 0. Initial conditions {t — 0) are found by a conservation of 
energy argument where the initial radius and infall velocity are given by 

ri2(t = 0) = i?i + i?2 



nit = 0) = ^/ — - M2$(Mi,ri2)] , (7) 
rrir 

where R2 is the maximum radius of cluster M2. 

In order to test the accuracy of this approach to the dynamics of cluster mergers, we 
compare the semi-analytic values of ri2 and Vq with the results of an N-body simulation of a 
two-cluster merger model (Berrington 2000). For this comparison, we assume that the two 
clusters have equal masses. This is the worst case, because when one cluster is smaller and 
less massive than the other, it is better described as a point source. In the N-body simulation, 
each cluster is approximated by = 80, 000 equal-mass particles, with 50% of the particles 
assigned to galaxies in each cluster and the remaining 50% assigned to a non-luminous dark 
matter which follows a King-model distribution with Wq = 6.25 (King 1966). The kinetic 
energy of the particles in each cluster is assumed to be in virial equilibrium with respect to 
the total gravitational potential energy prior to the interaction. Each cluster consists of 50 
galaxies with masses that are selected from a Schechter luminosity function with a power-law 
index of —1.25. The top panel of Fig. 1 shows the initial configuration of the merger event, 
and the lower panel gives the relative difference between the centers of mass calculated in 
the N-body simulation and semi-analytically. 
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The calculation ends at the collision time ^cou, which is defined as the time when the 
centers of mass of the two clusters coincide. As can be seen, the semi-analytic value of 
the Vo is accurate to within 10% at times preceding tcou- After tcou, the forward shock 
velocity, Vi, can no longer be approximated by equation 7 because the shock decouples from 
the gravitational infall of the merging cluster, and is assumed to be constant as the shock 
propagates outward from the cluster center. 

In this paper, we model the cluster merger shocks rather than the accretion shocks. An 
infalling cluster is assumed to be a virialized, quasi-spherical distribution of galaxies, ICM, 
and dark matter. A cluster merger shock forms when the cluster merger velocity exceeds 
the sound speed of the ICM. In contrast, an accretion shock is the result of a spherically 
symmetric infall of matter from the surrounding volume onto a cluster. In the context of 
the previous work by Fujita & Sarazin (2001), our definition of merger events includes not 
only their definition of a cluster merger event(M2 > O.6M1), but also their definition of semi- 
merger events (O.6M1 > M2 > O.lMi). In their calculations, accretion events can account 
for only ^ 10% of the nonthermal photon luminosity. Although each cluster of galaxies is 
at the center of an accretion flow, we do not consider nonthermal particle production from 
accretion shocks because their nonthermal photon production is negligible in comparison to 
cluster merger events. Note that supersonic cluster merger shocks are less likely to form 
at high redshift for two clusters of galaxies with given masses. This is because the relative 
velocities of the merging clusters will be associated with orbits with smaller separations at 
earlier times, and will therefore have smaller relative velocities. In reality, the mean cluster 
masses are smaller at early times. These lower mass clusters have a lower virial temperature. 
If this is considered in the calculation of the strengths of the shocks seen in cluster merger 
events, it is possible that stronger shocks may be seen at higher redshifts. 



Shock formation and nonthermal particle production depends on the properties of ICM, 
which are largely determined by the mass of the cluster and cosmological epoch. Below we 
describe the equations used to describe the properties of the cluster environment. 

X-ray observations of galaxy clusters show that these clusters are permeated by a hot 
ICM of temperatures in the range of 5 keV-12 keV. At these temperatures the gas is well 
described by an ideal gas. The sound speed of an ideal gas in terms of the gas temperature 
kTx is simply 



2.1.2. The Cluster Environment 
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where F represents the ratio of specific fieats wfiich we take to be equal to 5/3, is tfie mean 
atomic mass, and nip is the proton mass. 

The temperature of the ICM depends on both the mass of the cluster and redshift z. 
In order to calculate the cluster X-ray luminosity which gives the energy of X rays that can 
be Compton scattered, we first calculate the expected temperature of the ICM gas. We use 
the cluster M-T relationship from Bryan & Norman (1998), given by 

kTxiz) = 1.39/t, (i^) ' (^oh'E'{z)fs keV (9) 

where M is the mass of the cluster, h is the parameter defined hj H = 100 h km s^^ Mpc^^, 
the function E^{z) = n^il + zf + ^^{1 + z)'^ + and fr^ is a normalization constant 
taken to be 0.8. We have made use of the following standard definitions: 

with Ho as the current measured Hubble constant, pc is the critical density required to 
close the universe, A is the cosmological constant, and R is the current curvature radius of 
the universe. We force the constraint Qq + + Q\ = 1 as mandated by an infiationary 
cosmology. For all models in this paper, we set Qr = for a fiat universe. The function A^, 
is the critical density factor and marks the point at which a region of overdensity makes the 
transition from linear density growth behavior to a nonlinear density growth regime. This 
marks the birth of a cluster or the beginning of the initial collapse of a density perturbation 
out of the Hubble flow. Bryan & Norman (1998) calculate it to be 

ISn^ for Qo = 1 

Ac=-( 187r2 + 82x - 39x2 for = , (11) 
187r2 + 60x - 32^2 for = 

where x = Q(z) — 1, and Q{z) = flo{^ + zY/E'^[z). These relations are accurate to 1% in 
the range Q{z) = 0.1-1.0. 

A number of radiation fields provide seed photons for the Compton scattering by the 
nonthermal electrons present in the cluster. For the X-ray radiation field, we calculate the 
X-ray luminosity of a cluster of mass M by the observed relation (Arnaud & Evrard 1999) 

Lx{Tx) = 2.SSh-' X 10^^ (f^) ' ergs s'^ , (12) 

where kTx is defined by equation (9). Assuming the X-ray luminosity is uniformly distributed 
within a spherical cluster of radius R and an escape time of '^R/c, the mean energy density 
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of the X-ray photons is 

Ux(Tx) ^ 1.5/i^ X 10-^° 



Ri 



1 Mpc 



kT 



6 keV 



MeV cm 



-3 



(13) 



The thermal X-ray bremsstrahlung is approximated by a blackbody of temperature Tx in 
our calculations of Compton scattering. This is not a crucial approximation, because the 
X-ray energy density is ~ 2-3 orders of magnitude less than the CMB energy density, and the 
importance of Compton-scattered X-ray photons is furthermore reduced due to the Klein- 
Nishina effects. 

A second component that contributes to the mean photon energy density of the ICM is 
the stellar photon field. We assume that the galaxy luminosity function is well approximated 
by the Schechter luminosity distribution: 



ngal(iv) 



k 




—a 


L ' 




m 


exp 





(14) 



where n^axiL) dL is the differential luminosity function and is defined as the number of 
galaxies within the luminosity range L and L -\- dL. The parameter L* is the "characteristic 
luminosity" of a galaxy which marks the transition between the faint-end power law of slope 
—a and the high-end exponential cutoff. The constant normalizes the Schechter function 
to Ngai when integrated over all luminosities. We can use equation 14 to calculate the total 
stellar luminosity of 

/•oo 

L,tar = / dLL n^,,{L) = ^,T[2 - a]L* (15) 

where r[2 — a] is the gamma function. Typical values for L*, and a are given, e.g., in 
Paolillo, et al. (2001), and our adopted values are a — 1.1, a value of L* = 7 x 10^^ ergs s~^ 
corresponding to an absolute magnitude M* — —20.5 in Bj CoUess (1989), and = 100. 



Using the same approach as for the X-ray photon field, the estimated mean stellar 
radiation field becomes 



star 



1.4r X 10" 



1 Mpc 



MeV cm 



-3 



(16) 



The spectral distribution of elliptical galaxies is well approximated by a K3111 star (Pierce 
& Berrington 2002). Since higher density regions in rich galaxy clusters are dominated 
by elliptical galaxies, we will assume that the stellar radiation energy distribution can be 
approximated by a blackbody spectrum resembling that of a K3I11 star at temperature 
Tstax ~ 4000 K. 
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The energy density of the CMB is 

UcMB = 2.5 X 10-^(1 + z)^ MeV cm-^ (17) 

The CMB is the dominant contributor to total cluster photon energy density by approx- 
imately 3 orders of magnitude. The total cluster radiation field is simply the sum of the 
equations (13), (16), and (17). 



2.1.3. Cluster Redshift 

The cluster environment described in section 2.1.2 depends on the redshift of the cluster. 
The redshift evolution of the cluster depends on the chosen cosmological model. We have 
chosen to uniquely specify the cosmological model by the following parameters: critical mass 
fraction fl^, the curvature flu, the dark energy fl^, and the Hubble constant Hq. For this 
paper we will consider the model {flo, Qr, Oa) = (0.3, 0.0, 0.7). 

The observed redshift of a cluster is determined from the time since the Big Bang by 
the cosmological equation (Peebles 1993) 

^^^-Ho{l + z)E{z) (18) 

where E'^{z) is the dimensionless parameter defined following equation (9). Typical propa- 
gation times for a shock front to traverse a cluster are ~10^ yrs. This is a significant time 
span in that the evolution of the cluster environment due to cosmic expansion must be con- 
sidered. Cluster merger events can be initiated at any time throughout age of the universe, 
and evolving the cluster environment is necessary to accurately calculate the nonthermal 
photon and particle spectra. 



2.2. Nonthermal Particle Production 

Henceforth we denote particle kinetic energies by K^^p = me^p{'je^p — l)c^, total energies 
by Eg^p, and dimensionless total energies by e^, with the subscripts e and p referring to 
electrons and protons, respectively. First-order Fermi acceleration at a cluster merger shock 
is approximated by injecting power-law momentum spectra for the electron and protons. 
The total particle injection function in terms of kinetic energy is given by 

Qe,p{^e,p, t) = Ql^p[Ke,p{Ke,p + 2me,pC^)]~~^ {K^^p + nie,pC^) CXp 



K, 



e,p 



(19) 
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where s{t) is the injection index, and Xmax is the maximum particle energy determined by 
three conditions: the available time to accelerate to a given energy since the beginning of 
the merger event; the requirement that the particle Larmor radius is smaller than the size 
scale of the system; and the condition that the energy-gain rate through first-order Fermi 
acceleration is larger than the energy-loss rate due to synchrotron and Compton processes. 
The duration of particle injection is determined by the dynamics of the system, and the 
shock speeds are determined by equation (4). The shocks terminate at the time tacc,i for the 
forward shock and iacc,2 for the reverse shock when the shock propagates to the outer part 
of the cluster. The time t^cc represents the greater value of two acceleration times tacc,i and 
^acc,2- The constant normalizes the injected particle spectrum over the volume swept 
out by the shock front, and is determined by 



(iXe,p Ke,pQe,p{Ke,p, t) = vLmpVi{tf {nicuit)) , (20) 



where A is the area of the shock front, 1]^^^ = 1.2 is an enhancement factor due to the presence 
of ions heavier than Hydrogen, and (nicM(^)) is the number density of the gas averaged over 
the area of the shock front at time t. We assume a constant efficiency factor r]e,p = 5% 
for both electrons and protons, though r]e,p would depend on Mach number through the 
uncertain wave generation processes at the shock front. ^ The particle injection index, s{t), 
is determined by the shock properties described in §2.2.1. Typical values for the particle 
injection index are s ^ 2.5. 

The temporally evolving particle spectrum is determined by solving the Fokker-Planck 
equation in energy space for a spatially homogeneous ICM, given by 

dN{K,t) 1 92 , ^^.^^ d 



dt 2 dK"^ 



[D{K, t) N{K, t)] - — [[ku,t{K, t) - A,,,{K, t)) N{K, t) 



i=pp,p^,d ^ ' 

The quantity Xtot represents the total kinetic-energy loss rate, and Atot represents the total 
energy gain rate. Electrons experience energy losses from synchrotron, bremsstrahlung, 
Compton, and Coulomb processes. The total energy loss rate for protons is due to the 
effects of Coulomb processes. In addition, protons experience catastrophic energy losses due 
to proton-proton collisions [i = pp), proton-7 collisions (i = pj), and diffusive escape (i = d) 
on the time scale Ti{Kp,t). The spectra of the secondary electrons are calculated from the 



quantitative prescription for the shock efficiency as a function of the injection index is given by Miniati 
et al. (2001); see also Keshet et al. (2003) on this issue. 
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pion-decay products (see Appendix A. 3), and are subject to the same energy losses as the 
primary electrons. The calculation of the energy loss rates accounted for by the quantity X^ot 
for both electrons and protons is presented in detail in Appendix A. The nonthermal energy 
spectra resulting from the temporal evolution of equation (21) will produce a nonthermal 
photon spectra. The method used to calculate the synchrotron, bremsstrahlung, Compton, 
and TTo 7-ray radiation components to the photon energy spectra is presented in detail in 
Appendix B. 

The integration of the partial differential equation (21) poses a significant computational 
problem when considering the ranges of energies modeled in this paper. The problem stems 
from the diffusion that occurs at low energies, and the large energy losses that occur at high 
energies. We found that equation (21), when solved with a Crank-Nicholson implicit finite 
differencing method, is unstable to an oscillation instability at the high and low energies due 
to the large values that the convection derivative takes at both energy extremes. To solve 
equation (21), we settled on an implicit, upwind finite differencing method. This has the 
advantage of solving the oscillation instability problem, though with the disadvantage of a 
decrease in numerical accuracy for a given stepsize. 



2.2.1. Particle Source Function 

Nonthermal particle distributions are injected into the ICM according to the function 
given by equation (19). The normalization constant Q^p is calculated as described in §2.2. 
The particle distribution functions have spectral index s{t) and exponential cutoffs at kinetic 
energy Kj^ax, which are determined by shock acceleration theory as described below. 

Nonthermal particles gain energy through first-order Fermi shock acceleration at the 

rate 

K,„(t)^ 100^^ MeYs-' (22) 

(Sturner, et al. 1997; Reynolds 1996), where v^{t) = Vi{t)/{10^ cm s~^) is the shock velocity 
in units of 10^ km s~^, / is a factor giving the mean free path of a particle in units of its 
gyroradius, and the quantity rj accounts for the orientation of the shock front relative to 
the magnetic field. In our calculations, we set / = 10 and rj — 1. 

The first constraint that we impose on maximum particle energy is the available time to 
accelerate particles to a given energy. This constraint is determined through the expression 



Jo 



(23) 
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where vi{t) is found from equation (4), and B is assumed to be constant throughout the 
cluster volume. This constraint is the most important limitation on proton energy for most 
of the lifetime of the shock front. 

Competition with energy losses sets a second constraint on maximum particle energy. 
By equating the electron energy gain rate, equation (22), with the sum of the synchrotron 
and Compton energy losses that are the dominant electron energy-loss processes at high 
energies, we obtain a second energy constraint for electrons, given by 



= 2.8 X 10^ 



(1 + + 



MeV. (24) 



where B^j is the magnetic field in units of 10 ^ G. Considering synchrotron losses only, the 
corresponding maximum energy for protons is is 



K^s.,2isyn) = 2 X 10^^^3J= MeV. (25) 



For the parameters used in our calculations, proton energy losses due to Compton and 
synchrotron processes are smaller than photomeson and photopair losses with the CMB. 
These latter processes could in principle constrain the maximum proton energy (Berezinskii 
& Grigor'eva 1988), but are found to be negligible for the parameters used in the calculations. 
This is because the time scale to accelerate protons to ultra-high energies (> 10^*^ cV) 
requires ^ 3 x 10^^[frj/B_-jvl{t)]{Kp/10'^^eV) yr, which is longer than the Hubble time tn 
and much longer than the time scale over which the merger shocks persist in the merging 
clusters (~ O.I^h)- Although photomeson production does not pose a significant limitation 
to particle acceleration, it is still treated as a catastrophic loss in the calculation. 

In order for particle acceleration to occur, the Larmor radius of a particle must be 
smaller than the size scale of the system; otherwise, the particle escapes the system before 
any significant energy gain can occur. The third constraint imposed on maximum particle 
energy is therefore obtained by requiring that the particle Larmor radius be less than the 
size scale of the shock front. This constraint becomes, for both protons and electrons, 

.^,3^^.if..(^)MeV 

where A = 2i?2 and e is the electron charge. In all of our calculations, the limitations 
imposed by -ft^max,! and -ftrmax,2 are more restrictive for the electrons. However, for protons 
-f^max,3 is the Strongest limitation for the highest energy protons at the later stages of the 
shock evolution. The maximum energy Kmax{t) in equation (19) is the minimum of the 
maximum energies calculated in equations (23)-(26). Note that the size scale limitation on 
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maximum particle energy assumes that the coherence length of the magnetic field is larger 
than the size scale of the shock. 

The particle index s{t) is determined from the coUision and sound speeds of the system. 
The collision speed Vi{t) is obtained by solving equation (4) with the initial radius of the 
interaction given by the sum of the assumed radii of the clusters and the initial collision 
velocity given by equation (7). We assume that the subcluster is a point mass and the 
dominant cluster is described by the isothermal beta model density profile given in equation 
(3). 

The shock compression ratio is calculated from the equation 

where A4i(t) = Vi{t)/cs is the Mach number of the cluster merger shock speed, as computed 
from equations (4), (8), and (9). The ratio of specific heats is represented by T = 5/3 
appropriate to a nonrelativistic mono-atomic ideal gas. In Figure 2, the time evolution 
of the compression ratio of a shock front in a merging cluster with masses M — lO^^M© 
and m — lO^^M© is graphed. The radii of the two systems are assumed to be 1.5 Mpc 
and 0.75 Mpc, and the onset of the merger is at a redshift oi Zi — 0.1. Prom the shock 
compression ratio, we calculate the injected particle index s from 

The calculated shock speeds of cluster merger shocks range from ~1.5 to 3 times the ICM 
sound speed, with power-law indices ranging from ~ 5 to ~ 2.3. Figure 2 also shows the 
evolution of the typical power-law index as a function of time. 

In the merger event between two clusters of galaxies, as described in §2.1, two shock 
fronts develop. The shock front propagating into the dominant cluster, which we call the 
forward shock, has Mach number A4i{t). Similarly, the shock front propagating into the 
merging cluster has Mach number A42{t)- Writing down the compression ratio for the 
forward {Ci{t)) and reverse {C2{t)) shocks from equation (27), we get 

ui(t) 4 , ^ / s U2(t) 4 

where Ui{t) and U2{t) are the velocities of the forward and reverse shocks, respectively. In 
addition, v and v are the velocities of the shocked fluid for the forward and reverse shocks, 
respectively. Noting that v = voit) — v where vq is the infall velocity of the merging cluster. 
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and solving for the Mach number, we get 



and 



M2{t) 




9 4 
4{vo{t)-vf 



(30) 



where Ci is the sound speed in the dominant cluster, C2 is the sound speed in the merging 
cluster, and the Mach number is defined to be A^i,2 = ■^1,2/01,2- The values of ci and C2 are 
calculated from equation (8) with temperatures given by equation (9) for clusters of mass 
Ml and M2. The value of v is calculated by solving 



where Uiit) and n2{t) arc the ICM number densities in the dominant cluster and merging 
cluster at the positions of the forward and reverse shocks, respectively, and are calculated 
according to equation (3). We take the mean mass per particle /xi = /i2 = 0.6mp in the two 
clusters. Equation (31) was derived from the shock jump conditions by equating the energy 
densities of the shocked fiuids at the contact discontinuity. The compression ratios are then 
calculated from equation (27). These values are plotted in Figure 2. 



We present results of the numerical model described in §2 to simulate the temporal 
evolution of nonthermal particles accelerated by shocks formed in merging clusters of galaxies. 
Prom these particle distributions, we are able to investigate the expected nonthermal photon 
spectra for energies ranging from ~10~^ eV to ~10^^ cV that cover astrophysically interesting 
frequencies. At photon energies ^ 4 x 10^^ eV (~ 10 MHz), ground-based radio observations 
are subject to the ionospheric absorption cutoff. At photon energies > 10^^ eV (depending 
precisely on source redshift) the universe becomes opaque to 7 rays due to pair production 
attenuation on the diffuse infrared radiation fields. The electron-positron pairs produced 
by the ~ 10^^-~ 10^^ eV provide a secondary source of synchrotron and Compton-scattered 
radiation on the Mpc galaxy-cluster scale, but this cascade radiation is found to provide only 
a minor nonthermal synchrotron component to the nonthermal UV/soft X-ray emission, 
and is not included. To calculate detectabihty of clusters with ground-based Cherenkov 
telescopes, 77 absorption of the high-energy 7 rays must, however, be considered. 



Ill nijt) _ 1 + SMi^jt) / vojt) - V 
Ii2n2it) ~ l + 3M2^{t) V V 




(31) 



3. 



Results 
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The energy-loss time scales of nonthermal electrons for various parameters assumed 
in cluster environments are shown in Fig. 4. At low energies, Coulomb losses dominate. 
A maximum in the energy-loss time scale is found in the energy range from ^ 1 MeV to 
several hundred MeV. At higher energies, Compton and synchrotron losses dominate. At 
the highest energies, the Compton loss rate declines due to the onset of the Klein-Nishina 
reduction in the Compton cross section, and synchrotron losses become the dominant loss 
mechanism (see also Sarazin (1999)). The electron loss-rate from double Compton scattering 
is always much smaller than synchroton losses in our study, and can be safely neglected (see 
A.1). 

Figure 5 shows the nonthermal particle energy spectra at different times during the 
merger of a dominant cluster of mass Mi — IO^^Mq and radius Ri — 1.5 Mpc with a 
subcluster of mass M2 = lO^^M© and a radius R2 — 0.75 Mpc, respectively. The quantities 
Ml, M2, -Ri, and R2 are defined in §2.1.1. The assumed core radius Tc for the dominant cluster 
is 250 kpc. The ICM density of the dominant cluster is assumed to follow equation (3) with 
a central number density uq = 10^^ cm~^and /? = 0.75. These parameters are consistent 
with the Lx-Tx relation (12). A uniform magnetic field of strength l/xG was assumed 
throughout the cluster. The galaxy cluster environment and its evolution is calculated using 
the cosmological parameters {flo,flji,fl\) = (0.3,0.0,0.7). All models assume a Hubble 
constant Hq = 70 km s~^ Mpc~^. We also assume a Hydrogen to Helium number ratio of 10. 

Irrespective of the specific parameters of the simulation, we sec that at energies < 
1 MeV, the electron and proton spectra display low-energy cutoffs due to Coulomb diffusion 
and energy losses (the thermal particle distributions are not shown) . Electrons with energies 
greater than ~ 1 GeV lose a significant fraction of their energy on Gyr time scales. The 
dominant energy-loss mechanism for these electrons is Compton scattering, as Thomson 
losses on the CMB dominate synchrotron losses when B <^ 3(1 -|- //G (see Fig. 4). The 
dominant energy-loss mechanism for protons is secondary nuclear production, which operates 
on a time scale of about 2tH for an ambient density of 10~^ cm~^. Consequently, the high- 
energy cutoff in the proton distribution function results from the available time to accelerate 
particles and, at the longest times, competition of particle acceleration with catastrophic 
losses. Even though particles are injected with a range of power-law indices, most of the 
injection power occurs when the spectral index is hardest, yielding a proton spectrum that 
can be accurately described by a single power law (compare Fig. 2) with a number index of 
—2.5 for the forward shock, and ~ —2.2 for the reverse shock. 

Separate radiation components produced by nonthermal particles are shown in Figure 
6 for a shock formed in the merger of lO^^M© and lO^^M© clusters, with the interaction 
beginning at Zi — 0.3. The system has been evolved for 0.8 Gyr (z — 0.22) in panels (a) and 
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(b), and for 3.42 Gyr {z — 0) in panel (c). Panels (a) and (b) correspond to times before 
particle injection has stopped. The cluster magnetic field strength B = 1.0/iG in panels 
(a) and (c), and B = O.l/xG in panel (b). The radiation components of the primary and 
secondary electrons are shown by the short-dashed and dotted curves, respectively, and the vr'^ 
decay 7-rays are shown by the dotted curves. The curves are labeled by the various radiation 
processes, and the total spectral energy distributions are shown by the solid curves. The 
synchrotron and Compton radiations from primary electrons dominate the emission at times 
before primary injection stops, though 7r° decay 7 rays make a dominant contribution above 
100 TeV. At late times after particle injection ceases, the coohng primary electrons provide 
the dominant contribution at the lowest radio frequencies and at optical/UV frequencies, 
but secondary electrons can make the dominant contribution at 1.4 GHz. At 7-ray energies, 
primary bremsstrahlung and secondary vr'^-decay 7 rays make the dominant contribution at 
late times. Note that structure in some of the emission components is due either to knock-on 
electrons or to the separate effects of the stellar. X-ray, and CMB Compton energy losses on 
the electron energy spectra. As can be seen, the intensity of the synchrotron radio emission 
is very sensitive to the magnetic field strength. 

Figure 7 shows the temporal evolution of the total photon spectra radiated by a shock 
formed in the merger of lO^^M© and IO^^Mq clusters. The photon spectra are shown at 
various times for shocks that begin at redshifts Zi — 0.03, 0.05, 0.1, and 0.3. Each shock 
is evolved to the present time {z — 0.0), and the radiation spectrum at the present time is 
represented by the sohd curve. Typical lifetimes for the shocks are ^ 1 Gyr. As a result, 
only shock fronts with 2; > 0.1 show late-time evolution of the total photon spectra. 

At early times well before particle injection ceases (< 0.5 Gyr), the luminosity increases 
as the shock front accelerates particles and the particles accumulate. During this epoch, 
the injection spectrum is very soft and nonthermal bremsstrahlung and 7r° radiation can 
dominate the production of photons at high energies because of the relatively few high-energy 
electrons. The peaks in the spectra ai t < 0.3 Gyr spectra are, from low to high energies, 
due to Compton, bremsstrahlung, and 7r° radiation. The enhancement at low frequencies 
is due to nonthermal synchrotron radiation. As the front ages and the electron spectrum 
hardens (0.5-1 Gyr), the number of high-energy electrons increases and the Compton process 
strongly dominates the nonthermal emission at infrared energies and above. 

At times between tcoU and tacc (between ^ 0.9 and 1.1 Gyr in Fig. 7), the shock front 
propagates to the outer regions of the clusters, and the particle injection rate decreases. 
The number of high energy electrons begins to decrease as their previous numbers are no 
longer replenished. The consequence of this behavior is a softening in the synchrotron and 
high-energy Compton emission. This is seen clearly in the Zi — 0.1 and Zi — 0.3 panels. 
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As the Compton emission declines, the 7 rays from the p + p — > ttq — > 27 process become 
increasingly evident at times ^1.1 Gyr. 

The shock front terminates at ^acc after propagating through the cluster. At these late 
times, particle injection ceases, and the remaining population of nonthermal particles evolve 
according to the energy loss and diffusion processes described in §§A.l & A. 2, with no freshly 
accelerated electrons to repopulate the particle distribution. The system quickly becomes 
devoid of primary electrons, so that the bremsstrahlung and pion emission dominate at 
hard X-ray energies and above. The most extreme case is shown in the Zi = 0.3 panel of 
Figure 7. The spectrum at the present epoch shows three maxima produced by synchrotron 
radiation at kHz frequencies, Compton scattering near 100 eV, and combined nonthermal 
bremsstrahlung and TTo-decay emission near 100 MeV. The plateau in the radio emission at 
late times is formed by synchrotron emission from secondary electrons (1.4 GHz corresponds 
to 6 X 10-6 eV). 

Also shown in Figure 7 is the thermal bremsstrahlung emission from the ICM using the 
form for the emission spectrum given by Rybicki & Lightman (1979). The temperature and 
the luminosity of the thermal bremsstrahlung is computed from equations in §2.1.2. Thermal 
bremsstrahlung generally dominates the emission at 10 eV - 10 keV energies. It should 
be noted that the relations for temperature and luminosity of the thermal bremsstrahlung 
component do not take into consideration the merging process of the system. Computer 
simulations of expected X-ray temperatures and luminosities of merging clusters indicate an 
increase in the thermal bremsstrahlung luminosity and temperature (Ricker 1998; Schindler 
& Miiller 1993). Therefore the thermal bremsstrahlung curves presented in Figs. 7 and 8 
very likely underestimate the true X-ray luminosity and temperature of the gas during the 
merger event. At times ^ ^aco the relation between the initial temperature and luminosity 
will depend on mass ratio of the merger components, the dark-matter density profile of the 
clusters, and impact parameter of the merger event (Ricker & Sarazin 2001). 

Figure 8 shows a calculation for the same system as in Figure 7, except for a cluster 
with a magnetic field of 0.1 /iG. We consider these two values because of the uncertainty 
regarding the mean magnetic fields in clusters of galaxies. Despite Faraday rotation mea- 
surements indicating magnetic fields strengths B > 1.0/iG (Clarke, Kronberg, & Bohringer 
2001), a comparison of radio and Compton luminosities suggest field strengths B pa 0.1/xG 
(Fusco-Femiano et al. 1999). Possible explanations for this discrepancy may arise from an 
observational biasing of the Faraday rotation measurement of the magnetic fields-weaker 
magnetic fields mean brighter Compton emission for a fixed radio luminosity-and the possi- 
bility that the distributions of the magnetic field and nonthermal particles are anticorrelated 
due to the accumulation of electrons in weak magnetic field regions (Petrosian 2001). The 
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major change between the two figures is an increase in the synchrotron radio emission, which 
is roughly proportional to the square of the magnetic field. The system with the lower mag- 
netic field has a slightly brighter Compton component. 

Light curves at various observing frequencies are shown in Fig. 9 using the parameters 
for our standard case given by Fig. 7 with B — 1.0//G, and a cluster merger beginning 
at Zi — 0.3. The light curves of the nonthermal radiation exhibit a common behavior 
independent of frequency. At early times, the spectral power rises rapidly as the clusters 
merge. The peak emission occurs when the centers of mass of the two clusters pass at tcoib 
after which the emission exhibits a slow decay and approaches a plateau at times t > tacc 
when particle injection has stopped. The rate of decay of the emission increases with radio 
frequency due to the stronger cooling of the higher energy electrons, so that the decay is 
slowest at lower frequencies. Synchrotron emission from secondary electrons forms the late- 
time plateaus at radio energies. This behavior is also apparent for the hard X-ray emission, 
although it is formed by primary bremsstrahlung and both primary and secondary Compton 
radiation at late times. At 7-ray energies, the 7r°-decay emission forms a plateau of emission 
that dominates soon after t — ^acc- 

In Fig. 10, we show the spectral energy distributions and spectral components for a case 
where the efficiency rjp = 10% for protons and r)e = 1% for electrons. An enhancement of the 
number flux of nonthermal protons compared to electrons is expected in nonrelativistic shock 
acceleration because protons, unlike electrons, resonate with Alfven waves at all energies 
(Baring et al. 1999), and the Larmor radii of protons in the tail of the Maxwcllian distribution 
of the background thermal plasma is greater than that of the electrons, so that protons can 
more effectively scatter across the shock front. This theoretical expectation is borne out 
by the large proton-to-electron ratio in the cosmic rays. The enhanced proton contribution 
produces structure in the 7-ray spectrum even at i ~ icou when the emission is brightest. 
The 7r°-decay feature will be clearly evident in the spectrum of merging clusters of galaxies 
for such a large proton-to-electron efficiency ratio. 

Calculations of the hardest particle injection spectral index Smin formed in cluster merger 
shocks are shown in Fig. 11 as a function of the larger mass Mi of the two clusters, with 
the subcluster mass M2 = lO^^M©. To derive these results, we used equations in §2.1.1 to 
calculate the dynamics of the clusters, equations in §2.1.2 to calculate the sound speed, and 
equations (27) and (28) to calculate the spectral index from the Mach number. We also 
assume that the onset of the merger begins at redshift Zi = 0.1; softer injection indices are 
obtained for mergers at larger values of Zi because of the smaller maximum separations d 
at earlier times, given by equation (1). We calculate Smin for various values of Tc and 
The values (rc,/?) = (0.05,0.8), (0.05,0.45), (0.5,0.8), and (0.5,0.45) roughly correspond to 
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the extrema in the range of these parameters measured for 45 X-ray clusters observed with 
ROSAT (Wu & Xue 2000). Also shown are values of s^m for (r^,/?) = (0.179,0.619), which 
are the average values of these parameters for the 45 X-ray clusters, and (rc, (3) — (0.25, 0.75), 
which are the standard parameters used in the calculations. 

The uppermost curves with large Tc values display a maximum. At larger values of Mi, 
Smin hardens due to the stronger gravitational potential. Although Smin spans a wide range 
of values, the value of Smin corresponding to the average values of and j3 lies in the range 
2.1 < Sjniii < 2.4, and more realistically in the range 2.2 < Smin < 2.4 because clusters with 
masses 3 x 10^^ -IQ^'^Is/Iq arc very rare. Because Zi = 0.1, there is a constraint on the available 
time for the interaction to occur. At low values of Mi, the gravitational interaction is small 
so that the interaction takes a longer time to occur. The features in the curves between 10^^ 
and lO^^'^M© are a result of this effect. 

4. Discussion 

We have presented results of a computer simulation designed to calculate the temporal 
evolution of nonthermal particles accelerated by a shock front formed between two merging 
clusters of galaxies. We then calculated the expected total nonthermal photon spectra re- 
sulting from synchrotron, bremsstrahlung, and Compton radiation from both primary and 
secondary electrons, and TTo-decay 7 radiation from accelerated protons. Our results apply to 
the shocks formed in the merger of two clusters of galaxies, and not to the accretion shocks 
formed by spherical infall of matter during the process of structure formation. Theoretical 
studies suggest that accretion shocks should exist (Bertschinger 1985; Miniati et al. 2000), 
and that each cluster itself will be the center of its own accretion flow. Although X-ray 
and radio observations (Kassim et al. 2001; Markevitch & Vikhlinin 2001; Fusco-Femiano et 
al. 2003; Bagchi 2002; Valtchanov et al. 2002) of halos and relics provide strong evidence 
for the existence of merger shocks, there is also observational evidence for accretion shocks 
(EnBlin et al. 1998, 2001; Slee et al. 2001; EnBlin & Briiggen 2002). We have focused on 
merger shocks in this study. Previous work by Fujita & Sarazin (2001) would classify our 
simulations as accretion events, and they showed the luminosity of these accretion shocks 
to be ~ 10^^ ergs s~^ when Mi = 1 x lO^^M©, which is consistent with the results of our 
simulations. 

Beppo-SAX and Rossi X-Ray Timing Explorer (RXTE) observations show the presence 
of a hard X-ray (HXR) excess at energies > 20 keV from A1656 (Coma) (Fusco-Femiano 
et al. 1999; Rephaeli, Gruber, & Blanco 1999), A2256 (Fusco-Femiano et al. 2000), and 
A3667 (Fusco-Femiano et al. 2002a). Our simulations support a Compton-scattering CMB 
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interpretation of the HXR (Rephaeli 1979) and EUV (Hwang 1997; EnBlin & Biermann 1998; 
Lieu et al. 1999; Atoyan & Volk 2000) emission. As shown in figure 7 and 8, the nonthermal 
radiation from the merger shock during the post-colhsion evolution (> tcou) shows an excess 
above the thermal bremsstrahlung in both the HXR and EUV regions of the electromagnetic 
spectrum. For the standard parameters used, the calculated luminosity of the nonthermal 
component is Lnt ~ 10*^^ for the HXR between the 20-80 keV energies. For the Coma 
cluster, BeppoSAX observations indicate that the HXR flux is 2.2 x 10"^^ ergs cm~^ s~^ which 
corresponds to an observed luminosity of ~ 10^^ ergs s~^ in the 20-80 keV. This observed 
value is consistent with our simulations for the merger model between the times icou and 
^acc- Our model does match well with the calculated mass ratio of the Coma cluster merger 
model (CoUess & Dunn 1996). In A2256, the observed masses of the merging clusters are 
1.6 X IO^^Mq and 5.1 x IO^^Mq (Berrington, Lugger, & Cohn 2002). The observed HXR flux 
is 1.2 X 10"^^ ergs cm^^ s^^ which corresponds to a HXR luminosity ~ 10*^^ ergs s^^(Fusco- 
Femiano et al. 2000). By scaling our models to the observed masses of the merging clusters, 
then expected HXR luminosity is ~ 10^"^ ergs s^^ which is consistent with the observed values 
seen in A2256. 

The Extreme Ultraviolet Explorer has detected EUV emission in the ~60-250 eV band 
for a number of clusters. These clusters include Virgo, Coma, Fornax, A2199, A1795, and 
A4059 (Lieu et al. 1996a,b, 1999; Lieu, Bonamente, & Mittaz 1999; Berghofer, Bowyer, & 
Korpela 2000b; Bowyer, Lampton & Lieu 1996; Bowyer & Berghofer 1998; Mittaz, Lieu, & 
Lockman 1998; Bonamente, Lieu, & Mittaz 2001; Kaastra et al. 1999). Note, however, that 
no detection of EUV emission, other than that which could be attributed to thermal tail X- 
ray emission, is reported from analysis of EUVE observations for Fornax (Bowyer, Korpela, 
& Berghofer 2001), A2199 and A 1795 (Berghofer & Bowyer 2002), and A4059 (Berghofer, 
Bowyer, & Korpela 2000a). The measured EUV excess in the 60-250 cV band for Coma and 
A1795 is 2 X 10'''^ ergs s^^(Licu et al. 1996a) and 3 x 10'^^ ergs s~^, respectively (Bonamente, 
Lieu, & Mittaz 2001). The calculated EUV luminosity of the nonthermal component from our 
simulations is Ljvt ~ 10^^ for these EUV energies, and is consistent with both clusters. As 
mentioned before, our simulations are a good match to the Coma cluster merger scenario. 
The masses of the merging clusters in A1795 are uncertain (Oegerle & Hill 1994), and 
conflrming the results of our simulations with the observed luminosity of the EUV excess 
of A1795 is diflicult. An EUV excess of luminosity ~ 6 x 10"^^ ergs s~^ was detected from 
the cluster A2199 (Kaastra et al. 1999). However, A2199 lacks strong evidence for a current 
merger event. If any mergers occurred in the formation history of A2199, they occurred long 
enough in the past to allow the cluster to revirialize. It is possible that an adjacent cluster 
A2197 is beginning to interact with A2199 (Oegerle & Hill 2001), but seems unhkely that 
it is the origin of the EUV excess in the central 400 kpc of the cluster. Depending on the 
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masses of the merging clusters, the decreased luminosity of the EUV emission in A2199 is 
consistent with a cooling nonthermal particle population. Radio and X-ray observations of 
the central regions of A2199 indicate it is complex (Owen & Eilek 1998), and it is possible 
other mechanisms are producing the EUV excess emission. 

We treated in detail the merger of a lO^^M© cluster with a lO^^M© cluster having values 

of Tc and f3 in equation (3) approximately equal to the mean values found in a survey of 
ROSAT clusters of galaxies (Wu & Xue 2000). We found that the hardest spectral index 
reached by the merger shock in this system is Smin ~ 2.2 for Zi = 0.1. Because and (3 
arc not correlated with cluster temperature Tx and therefore cluster mass, as can easily be 
seen from the results of Wu & Xue (2000), our model results are reasonably representative 
of typical observed cluster mergers. 

We consider the detectability of nonthermal synchrotron emission with the Low-Frequency 
Array (LOFAR^). LOFAR is a multi-element interferometric low- frequency radio telescope 
designed for high resolution radio imaging in the 10-300 MHz frequency range. Figure 9 
includes several frequencies observable by LOFAR showing the optimization of the LOFAR 
telescope for detecting these objects. Estimated sensitivities for the LOFAR telescope range 
from 1.6 mJy at 15 MHz and 40 //Jy at 110 MHz. The quoted sensitivities are for one 
hour exposures for a single beam for a total of 13,365 dual polarization dipoles optimized in 
the 10-90 MHz, and 213,840 dual polarization dipoles optimized for the 110-220 MHz each 
with a maximum baseline of 400 km. With station diameters of 65m, resolution of LOFAR 
is approximately 12" at 15 MHz, and 2.4" at 110 MHz at the full extent of the array. Be- 
cause multi-element interferometric radio telescopes have trouble detecting large, exteneded 
sources, we must consider the angular size of the shock front. Given a size scale of ~1 Mpc, 
these features will easily be resolved at all frequencies visible to LOFAR at a distance of 100 
Mpc. As a worst case scenario, we will consider the extent of the array that is only capable 
of observing the shock front as a point object. For all but the most distant shock fronts and 
the lowest frequencies, central 2 km of the LOFAR array, known as the virtual core, is the 
optimum configuration. Within this virtual core is found 25% of the total collection area 
of the LOFAR array. Approximate peak luminosities taken from Figure 9 are 4 x 10'^ Jy 
at 15 MHz, and 5 x 10^ Jy at 110 MHz. For a shock front with an angular size of ~ 30' 
this gives ~ 1 mJy/n"and ~ 100 /xJy/D'' respectively. At these intensities, LOFAR will be 
able to detect any cluster merger shocks out to a distance of ~ 2000 Mpc in the 15 MHz 
frequency band and ~ 700 Mpc in the 110 MHz frequency band. These calculations neglect 
any absorption due to the interstellar medium within the Milky Way galaxy. 



^ http : / / www. lof ar .org/ 
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We use results derived using our standard system to estimate the detectability of cluster 
mergers with EGRET and GLAST.^ Generalizing the approach of Dermer & Schlickeiser 
(2002) to extended sources, the number of background counts Ni,g from extragalactic diffuse 
7-ray background photons with energies > detected within solid angle AO of a high- 
latitude source is given by 

O.B — do — i- \ -C/100 / 

where E^q = 100 MeV, Kb = 1.72(±0.08) x 10"^ photons cm-^ s-^ MeV-^ and as = 
2.10(±0.03) (Sreekumar et al. 1998). The on-axis effective area of the telescope is approxi- 
mated by the expression Ao(-E'-y/-E'ioo)"°, and At is the exposure-corrected time on source. 

The number of source counts produced by photons with energies > E^ and detected 
within solid angle AQ centered on the direction to a source is given by 
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(33) 



Here 6d{E) is the energy- dependent angular resolution (point spread function) of the tele- 
scope, 9s is the angular radius of the source, assumed to be energy-independent, and (f)s{E) is 
the source flux, which we write as 4>s{E) = (0o/i?ioo)(-E'/£'ioo)~"'' [photons cm^^ s^^ MeV~"^], 
with E in units of MeV. Writing the the uL^, spectral power as uL,y = Lq{E' / E[qq)"''' , where 
E' — {1 + z)E, E[qq = (1 + 2;) £^100, and o;,^ = 2 — is the uL^, source spectral index, we have 
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(34) 

The luminosity distance dL — lOOdioo Mpc, and the brightest luminosity reached at 100 MeV 
is Lq — 10^^L43 ergs s~^, noting from the calculations in Fig. 7 and 8 that with L43 1. 

The maximum sensitivity occurs when AQ = tt61{E) + tt9'^, giving 

[i-ao-a„) v-c/ioo/ 

Writing AO, — nO'^ — 9.6 x 10~^^2(°), where ^(°) is the telescope's acceptance angle in 
degrees, the background count rate 

;V,(> E,) . 1.65 X 10- '''"■^t"'' f^)"" ' ■ P'') 

i.i — ao V-c/ioo/ 



^ http : / /glast . gsf c . nasa. gov / 
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The sensitivity of detection at significance Ua is tfierefore given by 

. - .2.6X10-'^ ^'<-) (I^) . (37) 



^/Wg ' dj^o 0{°){1 - ao - au) \Eioo 

High confidence identification of a source requires that A^^ ^ 1 and n^j > 5. 

For the EGRET telescope on the Compton Gamma Ray Observatory^ Aq = 1200 cm^, 
Oo = 0, and a nominal observing time is At = lO^te s, with tg = 1. A source with a ra- 
dius of 0.5 Mpc will have an angular radius Og ~ 0°.3/(iioo- The point spread function of 
EGRET is 9d{E) = h°{E / Eio^f/'^ , which is larger than the source angular size except for 
the nearest sources or when observing at the highest energies, in which case the number 
of sources photons is small. Substituting into equations (36) and (37) gives Ns{> E^) ~ 
muU:i{E^/E^^oT'^-^/[{l-Oi.)dU] andn, ~ l.%^^U^{Ej E^o^Y'^+^-^'> /[{l-a,)dl^^]. Un- 
less Li3 ^ 1, we conclude that it is unlikely that EGRET could have detected nonthermal 
emission from merging clusters of galaxies. 

This is in disagreement with the estimate of Totani & Kitayama (2000), who argue that a 
large fraction of the isotropic unidentified EGRET sources could be due to the radiation from 

merging clusters of galaxies. Totani & Kitayama (2000) make very optimistic assumptions 
about the hardness of the spectral index of injected nonthermal electrons, and in fact claim 
that the merger shocks have such large Mach numbers that the injection index is nearly equal 
to 2.0. As shown in Fig. 11, such hard spectral indices are only possible in the most centrally- 
peaked clusters, but that the average cluster merger shock as inferred from observations of 
galaxy clusters has a minimum injection index between 2.2 and 2.3. Because the energies of 
electrons which Compton scatter the CMB radiation to > 100 MeV energies exceed ~ 200 
GeV, the steeper injection spectrum reduces the available power in these electrons by 1-2 
orders of magnitude, thus accounting for the difference between our expectations and the 
claims of Totani & Kitayama (2000), Kawasaki & Totani (2002), and Colafrancesco (2001). 
Although it is possible that a few of the unidentified EGRET sources are associated with 
merging clusters of galaxies, these would involve the less frequent events involving collisions 
of clusters with masses near IO^^Mq, and those with hard spectral indices or dark matter 
density profiles with strong central peaks. 

We have not taken into account efi'ects of a merger tree on the evolution of the nonther- 
mal particle distribution (Gabici & Blasi 2003). The protons that radiate in the EGRET 
and GLAST ranges will accumulate, and could contribute factors of 2-4 enhancements due 
to successive mergers, but not enough to significantly alter our conclusions. Our predictions 
are in accord with a recent analysis of EGRET data (Reimer et al. 2003) (see, however, 
Scharf & Mukherjee (2002)) that fails to detect gamma-ray emission from galaxy clusters 
that emit the most luminous thermal X-ray emission. 
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For the Large Area Detector on GLAST, Ao = 6200 crn^, ao = 0.16, and the on- 
source observing time for a scanning mode lasting tyr years is At = 0.2 x 3.16 x lO^tj^^ 
s, where the factor 0.2 is the exposure correction. The point spread function for GLAST 
is 9(1 = 3°.5{E / Eiqq)~'^/^ , which still generally dominates the source angular size except 
at the highest energies and for the closest sources. From this we obtain A^s(> -E^) ~ 
1300 V^43(^7/^ioo)"''-°-'7[(0-84 - a.)rf?oo] and ~ 14 L43x/V(^7/^ioo)"^+°-V[(0-84 - 
Q'f)'^ioo]- These estimates indicate that clusters of galaxies within 200 Mpc will be de- 
tected with GLAST, but we predict far fewer clusters of galaxies than predicted by Totani 
& Kitayama (2000) for the reasons discussed above. Note that the detection significance is 
only weakly energy-dependent, so that the criterion for detection becomes the requirement 
of detecting a few photons which is more easily satisfied at lower energies. 

Our results also have relevance to suggestions that clusters of galaxies contribute to 
the diffuse 7-ray background. Loeb & Waxman (2000) have argued that Compton-scattered 
CMB radiation from relativistic electrons accelerated by shocks formed during structure for- 
mation can make a significant contribution to the diffuse extragalactic 7-ray background, 
which has a featureless power-law spectrum between a few MeV to ~ 100 GeV (Sreekumar 
et al. 1998). If this emission is provided by cluster merger shocks, then our results indicate 
that this contribution cannot be significant for two reasons: first, the density distributions 
inferred from X-ray observations of thermal bremsstrahlung in Abell clusters imply an aver- 
age minimum injection spectral index that would produce Compton-scattered CMB emission 
that is softer than the spectrum of the diffuse 7-ray background. Although those few clusters 
with the hardest indices would contribute the most emission above 100 MeV, the superpo- 
sition of the emission from the softer sources would produce a concave spectrum at lower 
energies, which is not observed. Secondly, an associated hadronic signature would be ob- 
served even in the case that protons are accelerated with the same efficiency as the electrons 
at all times except near tcou (see Figs. 7 and 8). The hadronic feature would, however, be 
seen in the time-averaged 7-ray spectrum, so that calculations omitting hadronic emissions 
are incomplete (Keshet et al. 2003). If, as argued in §3, the hadronic acceleration efficiency 
exceeds that of electrons, then a 7r°-decay signature should be apparent in the background 
radiation (see Fig. 10) unless accretion shocks make the dominant source of high-energy 7 
ray cluster emission. In this case, most of the 7-ray emission should be located on the pe- 
riphery of the cluster, and the ttq 7-rays would be reduced dramatically due to the decreased 
lower gas density in the proximity of the protons (see Miniati (2003) for a spatially resolved 
model including both proton and electron emissions). 

Should it be the case that cluster merger shocks do produce a large fraction of the diffuse 
extragalactic 7-ray background, then both low hadronic efficiencies and dark matter halos 
with strong central density peaks in the majority of rich clusters would seem to be required 
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(see Fig. 11). Under the assumptions of isothermality assumed by the violent relaxation of 
dark matter halos which collapse from uniform expansion of the Hubble flow, dark matter 
halos should not have central density cusps. Furthermore, X-ray temperature and luminosity 
profiles arc inconsistent with strongly peaked central dark-matter density profiles (Makino, 
Sasaki, & Suto 1998; Wu & Chiueh 2000) despite the predictions of N-body simulations 
(Navarro, Frenk, & White 1997). Clusters of galaxies could still make a small contribution 
of order ~1-10% to the extragalactic diffuse 7-ray background radiation (Colafrancesco & 
Blasi 1998; Miniati 2002). 

Our simulations show that protons may reach ~ 10^^ eV for standard parameters, but 
are not accelerated to higher energies, in agreement with the gyroradius estimates of Norman, 
Melrose, & Achterberg (1995) for maximum cosmic ray energies. The maximum energies of 
protons are primarily constrained by the available acceleration times at early times. At 
later times, the proton Larmor radii may exceed the size scale of the shock, so that the 
protons then diffuse into the surrounding galaxy cluster of galaxies and undergo no further 
acceleration. Competition of the acceleration and p'-f energy-loss rates plays a small role in 
limiting the highest energy of protons for our standard parameters. 

For a diffusion coefficient corresponding to a Kolmogorov spectrum of turbulence, > 

10^^ eV protons can diffuse from the cluster on a Hubble time (Berczinsky, Blasi, & Ptuskin 
1997). Assuming Bohm diffusion, however, Berczinsky, Blasi, & Ptuskin (1997) and Volk, 
Aharonian, & Breitschwcrdt (1996) showed that clusters of galaxies provide a storage room 
for cosmic rays with energies < 10^^ eV on cosmological time scales. Our calculations 
confirm that ~10^^ eV protons remain within the cluster for times much greater than the 
lifetime of the shock front (see Fig. 5). For energies > 10^^ eV, the loss of protons from 
diffusion through the ICM will become significant over the lifetime of the shock front, and 
will become a restriction to the maximum energies. The maximum energy imposed by the 
Larmor radius size-scale criterion is always less than the maximum energy determined by 
the diffusion loss timescale for our models. For a given magnetic field, the Larmor radius 
of nonthermal particles is weakly dependent on the cluster mass (oc Af^/"'), and despite the 
possible range of masses for Mi, the maximum particle energy -Emax M^^^, and varies 
little over the mass range of Mi. Kang, Rachen, & Biermann (1997) argue that more rapid 
particle acceleration can take place in perpendicular shocks, but this requires preferential 
orientations of the magnetic field with respect to the shock direction. 
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5. Summary and Conclusions 

We have presented the resuhs of a computer code designed to calculate the nonthermal 
particle distributions of a shock front formed in the merger event of two clusters of galaxies. 
We have calculated nonthermal particle energy spectra for primary electrons and protons, 
and secondary electrons. Photon spectra were calculated for bremsstrahlung, Compton, and 
synchrotron processes as well as 7r°-decay 7 radiation from pp collisions. Our results apply to 
shocks that form at the interaction boundary between two merging clusters of galaxies, and 
not to cluster accretion shocks which form at the outer regions of the cluster. Evidence for 
nonthermal particle acceleration in merger shocks, which penetrate into the inner region of 
the accreting cluster, are provided by radio halos and relics and nonthermal X-ray emission 
from clusters of galaxies. 

The thermal X-ray bremsstrahlung was modeled from observations of luminosities, tem- 
peratures, and masses of clusters of galaxies. We modeled the nonthermal emission under the 
assumption that electrons and protons are accelerated with a 5% efficiency of the available 
gravitational energy that is dissipated during the course of gravitational interactions between 
two merging clusters of galaxies. Particle acceleration at shocks formed between merging 
clusters of galaxies produces a population of nonthermal electrons that Compton-scatters 
CMB photons. This process naturally produces a power-law distribution of photons in the 
40-80 keV energy range with luminosities consistent with HXR observations of the galaxy 
clusters A1656, A2256, A3667 and, possibly, A2199. In addition, the EUV emission observed 
in A1656 is also consistent with a nonthermal Compton-scattered CMB origin. Cluster mag- 
netic field strengths ~0.1/iG are in accord with this interpretation of the observed HXR 
excess and EUV emission. 

In contrast to the results of Totani & Kitayama (2000), Kawasaki & Totani (2002), and 
Colafrancesco (2001), we have argued that it is unlikely that more than a few of the isotropic 
unidentified EGRET sources can be attributed to radiation from nonthermal particles pro- 
duced by cluster merger shocks. Previous studies have assumed hard nonthermal spectra 
that give brighter 7-ray emission than implied by our simulation results. Such hard spectra 
can only be obtained in the infrequent merger events involving two very massive clusters, or 
between clusters where the dark-matter density profiles are centrally peaked. 

The diffuse extragalactic 7-ray background is a featureless power law with photon index 
of 2.10(±0.03). The dominant nonthermal 7 radiation components in cluster-merger shocks 
include Compton-scattered CMB radiation and bremsstrahlung from nonthermal electrons, 
and 7r°-decay emission from nuclear interactions involving nonthermal protons. The 7r° decay 
signature will be present if the hadronic acceleration efficiency exceeds the electron accel- 
eration efficiency, as expected in diffusive shock acceleration theory. Furthermore, spectra 
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calculated using parameters obtained from ROSAT observations of 45 Abell clusters (Wu 
& Xue 2000) produce spectra with slopes ~ 2.2-2.4. Unless dark matter density profiles 
are centrally peaked and hadronic acceleration efficiency is low, our results imply that non- 
thermal emission from merging clusters of galaxies can make only a minor contribution to 
the diffuse extragalactic 7-ray background. Using standard diffusive shock acceleration with 
mean ICM magnetic fields < l/^G, our results also indicate that the merger shocks in clusters 
of galaxies do not accelerate > 10^^ eV cosmic rays, though additional effects, such as shock 
obliquity and the presence of preexisting particle populations, could permit higher energy 
acceleration by these shocks. 

We wish to thank Armen Atoyan, Namir Kassim, and Paul Ray for useful discussions, 
and the referee for a detailed and constructive report. We would also hke to thank Steve 
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A. Particle Energy Loss Rates 

A.l. Electron Loss Rates 

Nonthermal electrons lose energy through Coulomb or ionization interactions with back- 
ground thermal electrons. Energy is lost through this process at the rate 

K^AE., t) ^ „,CM (^f^) [m-^), ( Al) 

where we have made use of the following definitions: 

S(t) = 27r-i / rfyy^exp(-|/), (A2) 

JQ 



This preprint was prepared with the A AS macros v5.0. 
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AW = 24-ln ,and (A3) 



,2 



The quantity njcM is the proton density of the ICM, /3eC is the electron velocity, Te{t) is the 
electron temperature given in Kelvin, and k is Boltzmann constant. Nonthermal electrons 
diffuse in energy space according to the diffusion coefficient for Coulomb scattering 

De,conl{Ke, t) = Ty^, TliCM ^^^^ ) kTe(t)E(t) (A5) 

(Miller, Larosa, & Moore 1996). 

For a fully ionized ICM, the bremsstrahlung energy loss rate (Blumenthal & Gould 
1970) is given by the expression 

/ gg6 \ 

KeMemiKe, t) = -T]^^ UiCM ( ) (1" 7e + 0.36) {K^ + m,C^) , (A6) 

where rj^^ = 1.3 is an enhancement factor for Helium that results from the Z{Z + 1) depen- 
dence of the loss rate. The synchrotron energy loss rate for electrons averaged over pitch 
angle in a magnetic field of strength B, Lorentz factor 7e, and velocity f3e is (Rybicki & 
Lightman 1979) 

Ke,.yn{Ee, t) = , (A7) 

where ot is the Thomson cross section. 

Electrons lose energy by Compton scattering ambient photons in the ICM environment, 
as described by equations (13), (16), and (17) for the cluster X-ray, stellar, and CMB radia- 
tion fields, respectively. The spectral distribution of each radiation field is approximated by 
a blackbody spectrum of temperature Tx (see equation (9)), 4000 K, and 2.7(1 ^ z) K at 
epoch z. The Compton energy loss rate is given by 

Ke„{Ee, t) = —car^lPl ^ UiaKf^{^Ji) (A8) 

i 

(Skibo 1993). The quantity (TRNile,^!) = <^c{le^i)/['^T{^ + le'&i)] is a Klein-Nishina cor- 
rection factor, where i?, = kTi/rrieC^ is the dimensionless temperature, and crc{€j) is the 
Compton scattering cross section for photons with energy = hu/nieC^ in the electron rest 
frame. 
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Electrons can also lose energy through the double (or "radiative" ) Compton process in 
which a second photon is created during an electron-photon interaction. Estimates by Ram 
& Wang (1971) and Gould (1975) show that the ratio of energy-loss rates through double 
Compton scattering and ordinary Compton scattering is ~l/4 when the typical photon 
energy in the electron rest frame is~ 50 MeV. This ratio only reaches unity when the rest 
frame photon energy is ~ 5 TeV. Synchrotron energy losses completely dominate at these 
energies, so the energy-loss rate through double Compton scattering can be safely neglected. 
A detailed treatment of electron energy losses in galaxy clusters is given by Sarazin (1999). 



A. 2. Proton Loss Rates 

The Coulomb energy loss rate for protons is given by 

Kp,cou\{Kp,t) = -7]^ nicM 



dx 



(A9) 



In addition to energy losses, the protons diffuse in energy space according to the diffusion 
coefficient 

(Huba 1994). 



Dp,coui{Kp, t) = ri" mcM -5 — ] kmt)E{t) (AlO) 

\ PpTTlpC 



~ 3 



Protons undergo proton-proton collisions to produce pions in which a proton loses 
of its energy per collision. After three collisions, the proton is assumed to be lost from the 
system. The timescale for this loss is 

Tpion{Ep,t) = — ^ . (All) 

C/3pniCM(o-H + 0.1(THe) 

(Sturner, et al. 1997), where ith and aue are the inelastic cross sections for pion-producing 
collisions with protons and helium, respectively. Measured values for these cross sections are 
an — 28 mbarn and ane — 100 mbarn well above threshold, respectively (Meyer 1972). 

High-energy protons interact with the CMB photons due to photomeson production. It 
requires about 5 interactions for protons to lose a substantial fraction of their energy due to 
interactions with photons with proton rest-frame energies in the range 200 MeV-500 MeV, 
and about 2 interactions for photons whose energies in the rest frame of the protons are 
greater than 500 MeV. The time scale for the catastrophic loss is given by 

c /"°° n (e) 

T^lAEp, ^) = 2 X (ri^ J_, ~ kp,{e')ap,^^{e') (A12) 
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where (Tp^-^T^{e') is the cross section of the proton-photon interaction, and e' = 7pe(l — A*) is 
the photon energy in the rest frame of the proton. Here /x = cos 9, and 9 is the angle between 
the directions of the interacting photon and proton. We let (Tp^_+7r(e') = 380 /ibarns and 120 
/ibarns, and the inelasticity coefficient kp^ = 0.2 and 0.6, for 200 MeV < e' < 500 MeV, and 
e' > 500 MeV, respectively (Atoyan & Dermer 2001). The number density of CMB photons 
at redshift z is denoted by n^(e). 

Protons diffuse through the ICM at a rate that depends on the the diffusion coefficient. 
This diffusion process is a random walk which allows particles to escape from the system. 
Protons that escape from the cluster merger are treated as a catastrophic loss with time 
scale (Volk, Aharonian, & Breitschwerdt 1996) 

TdiEp) ^ i??/KBohm(^p), (A13) 

where we employ the Bohm diffusion coefficient K.Bohm.{Ep) = (3cri,/3 for particle diffusion, 
and tl is the particle Larmor radius. 

A. 3. Pion Production from pp-CoUisions 

Nonthermal protons collide with ambient gas particles to produce secondaries through 
several channels. The five dominant channels are (1) p + p ^ 7r*^ + X;(2) p + p ^ 7r"'" + X;(3) 
p + p ^ TT" + X;(4) p + p ^ ir+ + A:;(5) p + p ^ i^" + X. Here "AT" refers to any other 
byproducts produced in the reaction other than the particle indicated. It is also understood 
that channel (2) does not include the contribution resulting from the production of deuterium 
in the reaction p + p — > tt"^ + d, which can be treated separately but is a small contribution 
to the total secondary production. The secondary electron production is dominated by the 
first two pion channels, and are the only two channels considered in our calculations. Our 
calculation technique follows the method described by Dermer (1986a,b) and Moskalenko 
& Strong (1998); sec also Blasi & Colafrancesco (1999). This model combines the isobaric 
model (Stcckcr 1970) at energies <3 GcV and the scaling model (Badhwar et al. 1977) for 
energies > 7 GeV. We use a linear combination to join the two models in the transition 
region between 3 GeV and 7 GeV. 

For proton energies <3 GeV, the production of pions produced by proton-proton col- 
lisions is mediated by the excitation of a A3/2 isobar. Assuming the outgoing A3/2-isobar 
of mass mA is coUinear with the initial direction of the colliding protons in the center-of- 
momentum system (CM) with isotropically distributed decay products in the laboratory 
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system (LS), the pion distribution in the LS of the isobaric model is given by 

(A14) 

where the Hcavyside function H[x; a, h\ is defined to be = 1 if a < x < 6; otherwise = 0. The 
function N^(=f) ~ I'k.l'-ni^-^ ^ ^"^^ is the pion LS kinetic energy. The forward (+) and 
backward (-) moving isobars Lorentz factors are = 7c7^(l ± /3c/3^), where 7c = y^/2mp 
is the Lorentz factor of the CM with respect to the LS, and 7^ = (s + + m^) /2s^^'^m/^ 
is the Lorentz factor of the isobar in the CM. The pion Lorentz factor in the rest frame of 
the A isobar is 7^ = {rn\ + — rrip) / 2171^111 

The pion spectrum resulting from the pp collisions is calculated by evaluating the fol- 
lowing integral over the isobar mass spectrum: 



Wr{Kp) / dmA B-^{m/\)f{KTr;Kp,mA)- (A15) 



By^{m/^) is the normalized isobar mass spectrum given by the Breit-Wigner distribution 

{rriA - {mA)) + 
with average isobar mass (mA), and normalization factor 

-I -1 



Wr{Kp) 



^fs-mp- {niA) \ . fnip + m^ - {mA) 

arctan ^ — ^ — arctan ' 



(A17) 



Our calculations only consider the A3/2(1232) isobar which has a resonance width F = 
1/2 X 115 McV (Particle Data Group 1984). 

For protons whose energies are greater than 7 GcV, we adopt the scaling model of 
Stephens & Badhwar (1981). The Lorentz invariant cross section given by Badhwar et al. 
(1977) and Stephens & Badhwar (1981) for neutral and charged pion production is given by 

^.-^ =^e?.(^p)(l-C)®exp[--Bp^/(l + 4mJ/s)], (A18) 

where we have made use of the following definitions: 

6;^±(Ep) = (l+4mJ/s)^ (A19) 
GAE,) = (1 + 22,Ef-'){l - ^mllsf, (A20) 

Q = (Ci - + C3pi)/yi + 4m2/s, (A21) 
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C.= sJxl + {Als){pl + ml), (A22) 

" s/ {s — — m\Y — 4ml.mx 

The constants A, B, Ci^2,3, and 71 are given by Badhwar et al. (1977) and Stephens & 
Badhwar (1981), and 9 is the pion LS polar angle. Adopted values for the constants A, 
B, Ci^2,3, and TZ are given in Table 1. The quantity mx depends on the channel under 
consideration and are (1) mx = 2mp, (2) nix = ''^p + nT-n, (3) mx = '^rUp + m-;^, (4) 
mx = mp + mn, and (5) mx = '2mp + m^- 

To calculate the LS energy distribution of pions, we integrate the following over the 
pion LS polar angle: 

where — 1 < cos6'max < 1 and is defined as 

cos 9^^ = (^^E^ - '-ff"" ) . (A25) 

The quantity {rjaj^^Ep)) is the inclusive cross section and its values are taken from equation 
(5-9) of Dermer (1986a). 

The pions will decay to produce muons. These muon are created fully polarized. This 
results in a decay asymmetry. The muons created by the pion decay are created in the 
rest frame of the pion with Lorentz factor 7^ = {m'^+m'j^)/2mT^m^ ~ 1.039, and (3^ ^ 0.2714. 
In the rest frame of the muons, the resulting electron decay spectrum is given by Commins 
(1973); Orth & Buffington (1976) 



dE* dcos9* m^ 



l+«l^|cos<)' 



(A26) 



where E* is the electron energy in the muon rest frame and helps to define the quantity 
e = 2E*/m^. The angle between the polarization angle of the parent muon and the resulting 
electron momentum is 9*. The quantity ^ = ±1 for e^. 

For a given electron energy E^, the normalized electron decay energy distribution in the 
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frame where the pions are isotropic is (Dermer 1986a) 



-(/'(<;_) - iIj{uj_) + (p{uj-) - 4>{',+), <;+ < a;_ < <?_ < cu+, 
t/'(a;+) - + - 0(<?+), <?+ < a;_ < a;+ < 

. 0, uj+< 



(A27) 



where we have introduced the following definitions: 



u± = 
(l){x) = 



7ji77rmji(l±/3,r) 

(1±/3m)/2 



.5 ^2 



87;la; 



(3-n/32)(l-/32) 



4x(3+/3g-4n/3g) 



1 



(A28) 



-U 9/^ In T - 6(m+2m/3^+3)x2 , 16(M+3M/3^+2)a; 
\^U-r«,;mx (l+/3^)2 3^4-/9,, 13 



where -u = and ,^ = 1 for positrons, and ,^ = — 1 for electrons. Figure 3 shows 

the energy distributions of electrons and positrons formed through pp collisions by mono- 
energetic protons with a variety of energies for channels (1) and (2) (compare Murphy, 
Dermer, & Ramaty (1987)). 

The resulting electron or positron spectrum formed by the decay of isotropically pro- 
duced pions is 

dN{E,-y,) 



Qe{Ee) 



dE, 



(A29) 



with g^(7^) given by equation (A24), and % = {EjE"^^ + E'^^''/E^)/2 if E^ > E^^ = 
m^(l 69.8 MeV; otherwise, = 1. 



A. 4. Knock-On Electrons 

Nonthermal protons interact with thermal electrons through Coulomb collisions to pro- 
duce a population of nonthermal knock-on electrons. The knock-on electron source function 
in electrons cm~^ sec""^ is given by 

Qkile) = 1.75 nicMir) 47r / d-fp %^e{le, Ip) -T^ (A30) 

J71 d-fp 

(Abraham, Brunstein, & Cline 1966) where the factor 1.75 accounts for the presence of 
heavier nuclei in the ICM and cosmic rays, 7e is the Lorentz factor of the electron, 7^ is 
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the Lorentz factor of the proton, and dJp/ d'jp is the differential cosmic-ray proton intensity 
in particles cm~^ sec~^ sr~^. The differential probability for the production of an electron 
having a total energy per unit rest-mass energy in d'ye at 7e by the collision of a cosmic-ray 
proton of Lorentz factor 7p in units of cm^ s~^ is given by 



^p^eile,lp) die 



IP 



(7e - 1)2 mp(7e - 1)72 



+ 



2771272 



d7e (A31) 



(Abraham, Brunstein, & Cline 1966) where r^. is the classical radius of the electron. The 
maximum transferable energy is given by 



7n 



= 1 + 



72 - 1 



(A32) 



The lower bound of the integration in equation (A30) is determined by solving the inequality 
7e < 7max for 7p. The result is 



> ^1 = 



7p ^ 71 



2m, 



■(7e - 1) + 



mz 



Ami 



If 



1/2 



(ASS) 



B. Photon Production 

Nonthermal photons are produced by electron synchrotron radiation, electron-electron 
and electron- nucleon bremsstrahlung, and Compton-scattered CMB, stellar, and X-ray ra- 
diation fields. These components are calculated for both primary electrons and secondary 
electrons and positrons. Neutral pion-decay 7 rays are also produced in proton-nucleon 
interactions. We discuss each component in the order mentioned above. 



B.l. Synchrotron Radiation 

We define the synchrotron radiation emissivity as Qsynch = dNgyj^ch/dV dt dE^. It is 
calculated from the expression 

\ de sin^e) / dKe JeiKe,t) F{u/iy,), (Bl) 

n umec J Jo Jk,,^,^ 

where 9 is the electron pitch angle, — hv is the photon energy of frequency i/c — 
0A2^'^B_7sm{9) in units of Hz, and Je{Ee,t) is the electron intensity (units of electrons 
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cm ^ s sr ^ E^^). The function F{v/vc) is defined through the expression 



F{v/vc) 



V 



dx K^i^{x), 



(B2) 



where K^i^ is the modified Bessel function of order 5/3. The integration is performed nu- 
merically with approximate expressions for F{u/uc) given by Ginzburg & Syrovatskii (1965). 
Electrons are assumed to be isotropically distributed in pitch angle. In addition, we in- 
clude the effects of free-free absorption of the synchrotron spectra in the source using the 
absorption coefficient (Ginzburg & Syrovatskii 1965) 



Q;ff(cm ^) = 10 



n 



ICM 



3/2 



17.7 + In 



n3/2' 



(B3) 



where the electron temperature Tg is given in Kelvin, the ICM particle number density is in 
units of cm~^, and the frequency u is in Hertz. 



B.2. Bremsstrahlung Radiation 

Wc calculate the electron-electron (e-e) and electron-nucleon {e-n) bremsstrahlung emis- 
sivity by the expression 

Qhrem,e-e,e-n{E-y,t) = 47rniCMA^^' / dK^ Je{Ke,t) —- , (B4) 

"'JSTe.mm / e-e,p 

where we have included the correction factor A^" for the presence of Helium. Approximate 
values for the Hehum correction factor are A^g = 1.2 and Agg = 1.4. The expressions for 
the differential cross sections are adopted from equation (Al) in Haug (1975) for electron- 
electron interactions and from equation (3BN) in Koch & Motz (1959) for electron-nucleon 
interactions. 



B.3. Compton Scattering 



We calculate the Compton scattering emissivity using the full Klein-Nishina cross section 
(Jones 1968; Blumenthal & Gould 1970) for relativistic electrons (7 » 1) by the expression 

nj{E,) 



Qcorap,j{Ej, i) 



STTrlirriecY 



dEt 



E, 



2qln{q) + {1 + 2q){l - q) + 



dE, E-^J^{E^,t) 

{qee,tni-q) ' 
2(l + ?ee,t) . 



(B5) 
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We define 



^ed{Ee — E^) 



1 + W1 + 



[rrieC 



2\2 



E^Et 



nj{Et 



exp 



- 1 



and 



(B6) 
(B7) 

(B8) 

(B9) 



The energy Et is target photon energy, and Ef'^^^{Et) is the minimum energy required by 
an electron to scatter a photon from Et to energy E^. The temperature, number density, 
and energy density of the j*'* photon component is given by Tj, rij, and Uj, respectively, 
from §2.1.2. 



B.4. Pion-Decay 7 Rays 

We use the method developed by Dermer (1986b) to calculate the 7r° emissivity from 
proton-proton coUisions. The 7r° emissivity is calculated by 

Qy{K^o,t) = 47rnicM H dK, J,{K,, t) (BIO) 

J thresh ClJX^O 

where we have defined the neutral pion kinetic energy as K^o, and the differential cross 
section for 7r° production is calculated from the cross section given by Dermer (1986b) which 
uses the isobar model of Stecker (1973) for proton energies less than 3 GeV and the scaling 
model of Stephens & Badhwar (1981) for protons of energies greater than 7 GeV. A linear 
combination of the two models is assumed between 3 and 7 GeV. The 7-ray emissivity from 
TTo decay is calculated from the expression 

Q17,{E,, t)^2r dE^. ^^y^^'f, (Bll) 

Considering collisions with heavy nuclei, Dermer (1986b) has shown that the 7-ray emissivity 
is enhanced by the multiphcative factor 1.45. The total 7-ray emissivity therefore becomes 

QX = l-45gy^ (B12) 



X [Code Units] 




Fig. 1. — Comparison between the semi- analytic model used in subsequent calculations and 
results of an N-body simulation. The initial configuration of an N-body simulation of a 
cluster merger between two equal mass clusters is shown in the top panel. The relative 
differences Ar and Av between the semi-analytic and numerical calculations of the center- 
of-mass distance and the infall velocity of the two clusters are shown in the lower panel. 




Fig. 2. — The compression ratio of the forward (dashed curve) and reverse (sohd curve) 
shocks as calculated by the method described in §2.2.1. The particle injection spectral index 
for the forward (dotted curve) and reverse (dot-dashed curve) shocks is associated with the 
compression ratio, and is calculated for a merger system consisting of a dominant cluster 
of mass IO^^Mq and a poor subcluster of mass lO^^M©. Initial onset of the shock begins 
at a radius of 1.5 Mpc from the cluster center with an initial redshift of Zi — 0.1. Particle 
acceleration occurs over a period t^cc = 1-1 Gyr for these parameters, after which injection 
stops. Acceleration at the reverse shock ends after the reverse shock passes through the 
merging cluster. 
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Fig. 3. — Energy distributions of secondary positrons {left panel) and electrons (right panel) 
produced in p-p nuclear collisions by protons with energies i^p(eV) = {3.16 x 10^,5.62 x 
10®, 3.16xl0^ 2.15x10^°, 1.0xl0l^2.5xl0l^6.3xl0l^l.58xl0l^ 4.0x10^1, 1.0x1024} (from 
bottom to top) for the p+p +X channel and protons of energies Kp = {1.0 x 10^, 2.15 x 
10^ 5.0 X 10^ 5.0 X 10^°, 1.0 X 10^2^ 2.5 x 10^^, 6.3 x 10^^ 1.58 x 10^^ 4.0 x lO^^, 1.0 x 10^^ eV} 
(from bottom to top ) for the p + p + X channel. 
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Ee [log (eV)] E, [log (eV)] 

Fig. 4. — The instantaneous energy-loss timescale for electrons in a merging cluster scenario 
at redshift z = 0.0. In the left panel, the proton density n = 10~^ cm~^ and B = 1/iG. 
Only the dominant CMB radiation field is considered. In the right panel, n = 10~^ cm~^ 
and B = 0.1//G. 
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Fig. 5. — Primary electron {left panel) and proton {right panel) energy spectra for a shock 
formed by a merging cluster using a standard parameter set with Mi = IO^'^^Mq and M2 = 
1O^^M0 and B = 1.0 /xG. The interaction begins at redshift Zi = 0.3 and is evolved to the 
present time {z = 0.0; solid curve) with spectra shown at the various times given in the 
legend. The redshifts corresponding to the time delays of 0.1, 0.5, 0.8, 1.0, 1.2, 1.5, and 3.42 
Gyr are z = 0.289, 0.248, 0.218, 0.199, 0.180, 0.153 and 0.0, respectively. 
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Fig. 6. — Separate radiation components produced by nonthermal particles accelerated by a 
shock initiated at Zi = 0.3 resulting from the merger between lO^^'Mo and W^Mq clusters. 
The separate radiation components are labeled in the figures. Emission from the primary 
electrons are represented by the dashed curves, and emission from secondary electrons and 
positrons and tt'^ 7-rays are represented by the dotted curves. The sohd curves are the total 
nonthermal photon spectra. Panels (a) and (b) show the nonthermal spectrum at 0.8 Gyr 
{z — 0.22), and panel (c) at 3.42 Gyr [z — 0.0) after the merger shock forms. Magnetic 
field strengths are B — 1.0//G for panels (a) and (c), and B — 0.1//G for panel (b). Panels 
(a) and (b) correspond to times before the termination of the nonthermal particle injection. 
Panel (c) is approximately 2.3 Gyr after particle injection has ceased. 
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Fig. 7. — Temporal evolution of the total spectral energy distribution for a cluster merger 
shock which begins at redshift Zi = 0.03 {top left), Zi = 0.05 {top right), Zi = 0.1 {bottom left), 
or Zi = 0.3 {bottom right). Each shock is evolved to the present {solid line). Intermediate 
times are indicated in each panel. In all cases the shock survives for ~ 1.1 Gyr. Cluster 
merger shocks that develop at redshifts > 0.1 do not persist to the present epoch. The 
thin solid curve shows the expected thermal bremsstrahlung for a IO^^Mq cluster, whose 
temperature is given by equation (9). The figures with two red curves represent the thermal 
bremsstrahlung at Zi (lower curve) and z — 0.0 (higher curve). The integrated luminosity of 
the thermal bremsstrahlung is 10^^ ergs s~^. 
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Fig. 8. — Same as Figure 7, except with an assumed cluster magnetic field of 0.1/i.G. 
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Time [Gyr] Time [Gyr] 

Fig. 9. — Light curves at various observing frequencies produced by a shock formed in a 
merger between 10^^ M© and IO^^Mq clusters that begins at Zi — 0.3 and is evolved to the 
present epoch {t — 3.42 Gyr). All light curves are for a magnetic field strength oi B — 1.0//G 
unless otherwise noted. Radio light curves in Jansky units are given at 15 MHz, 30 MHz, 
74 MHz, 110 MHz, and 1.4 GHz on the left panel, and light curves in energy flux units are 
given at 40-80 keV, > 100 MeV, > 1 GeV and > 100 GeV in the right panel. The 15 MHz 
light curve is also calculated with a magnetic field strength of 5 = O.l/xG {dotted curve). 
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Fig. 10. — Same as the Zi — 0.3 case in Fig. 7, except that r]e — 1% and r]p — 10%. Right 
panel shows the separate radiation components at time t — 0.8 Gyr using the same notation 
as Fig. 6, and the left panel shows the total spectral power at various times. 
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Fig. 11. — Calculations of the hardest particle injection spectral indices Smin formed in cluster 
merger shocks as a function of the larger mass Mi of the two clusters, for various values of 
Tc and f3. The values of Tc are given in Mpc. The minimum spectral index for the forward 
and reverse shock is shown in the figure. 
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Table 1. Adopted values for constants defining the invariant cross section for p + p — > X 

reactions 



Constant 


(1) 


Channel 
(2) 


(3) 


A 


140 


153 


127 


B 


5.43 


5.55 


5.3 


Ci 


6.1 


5.3667 


7.0334 


C2 


-3.3 


-3.5 


-4.5 


C3 


0.6 


0.8334 


1.667 


7^ 


2 


1 


3 



